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ELECTROMAGNETIC  FIELDS  IN  THE  OCEAN  NEAR  A  SHORELINE 


ABSTRACT 


The  electromagnetic  fields  in  the  ocean  resulting  from  ground  wave  excitation 
over  land  by  natural  electromagnetic  noise  such  as  Micropulsations  or  ELF  Atmos¬ 
pherics  are  investigated  analytically.  A  simplified  two  dimensional  model  is  used. 
The  purpose  of  the  study  is  to  determine  the  relative  importance  of  electromagnetic 
energy  propagation  via  air  into  the  ocean  versus  propagation  through  the  soil  and 
ocean  bottom  into  the  sea.  Vertical  profiles  of  the  horizontal  electric  field  are 
plotted  for  several  values  of  frequency,  ocean  depth,  and  the  distance  from  the 
shoreline.  The  relative  importance  of  the  soil  path  Is  seen  to  be  greater  for 
shorter  distances  from  the  shore  and  for  deeper  oceans. 


1. 


Introduction. 


The  purpose  of  this  report  is  to  present  the  results  of  an  analytical  study  on 
the  relative  importance  of  electromagnetic  energy  propagation  via  air  into  the  ocean 
versus  propagation  through  the  soil  and  ocean  floor  into  the  sea  employing  the  highly 
simplified,  two  dimensional,  model  shown  in  figure  1.  More  complicated  analytical 
models  which  may  represent  a  closer  approximation  to  the  real  three  dimensional  situ¬ 
ation  will  be  discussed  in  future  reports. 


In  the  present  model  all  fields  are  considered  to  be  independent  of  z(i.e. 
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and  the  magnetic  field  H  has  a  z  component  only.  The  complete  electromagnetic  field 

can  thus  be  expressed  in  terms  of  H  as 

a)  H  «  h 
z 


b)  E 


:)  E 


+  1  9h 


o-iu)c  3y 
-1  3h 


(1) 


y  cr-iwe  dx 

where  a  is  the  conductivity  and  e  is  the  permittivity.  The  permeability,  uq,  of 
all  media  is  taken  to  be  that  of  free  space.  A  time  factor  of  exp  (-iuot)  is  assumed 
throughout.  The  electromagnetic  noise  from  natural  sources  over  land  is  assumed  to 
produce  a  nearly  vertically  polarized  plane  wave  propagating  in  the  +  x  direction. 
The  field  established  in  the  ocean  when  this  wave  is  incident  on  the  shoreline  is 
to  be  determined. 

(1) 


The  incident  field  is  given  by 
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(3) 


In  order  to  obtain  the  field  in  the  ocean  vre  will  treat  the  problem  in  the 
following  manner. 

We  consider  the  field  to  be  a  superposition  of  two  fields,  one  being  the  field  due 
to  the  energy  entering  the  ocean  through  the  air-sea  boundary  and  the  other  Being  the 
field  due  to  the  energy  entering  through  the  horizontal  land-sea  boundary.  The  energy 
entering  through  the  vertical  land-sea  boundary  is  neglected  in  this  treatment. 

The  magnetic  field  on  the  surfaces  of  the  ocean  (both  top  and  bottom)  will  first 
be  obtained  under  the  assumption  that  the  ocean  is  a  perfect  conductor. 

To  find  the  field  at  some  point  in  the  ocean  the  assumption  is  made  that  wave  enters 
the  ocean  normal  to  the  surface. 


Contribution  through  Air-Sea  Boundary. 

In  this  section  we  shall  determine  an  approximate  expression  for  the  magnetic  field 
in  the  ocean  due  to  the  energy  entering  the  air-sea  boundary  only.  We  begin  by  deter¬ 
mining  the  magnetic  field  at  the  top  surface  of  the  ocean. 

Consider  equations  (2)  and  (4)  which  describe  the  incident  wave  in  the  air  above 
the  ground.  For  the  range  of  frequencies  of  interest  here  the  quantity  CU)eQ/o2]  «  1. 

We  note  then  that  the  magnetic  field  in  air  is  not  very  sensitive  to  a  change  in  con¬ 
ductivity  and  therefore  will  propagate  past  the  shore  line  relatively  unperturbed. 

< 
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3. 


Neglecting  the  attenuation  (i.e.  considering  a  m  ®)  we  have  the  following  approximate 
expression  for  the  magnetic  field  at  the  top  surface  of  the  ocean 
ik  x 

hts  *  e  (5) 

U) 

where  k  =»  — 

o  c 


The  assumption  of  infinite  conductivity  used  in  order  to  obtain  h  is  now 
dropped.  The  conductivity,  however,  is  still  considered  high  enough  so  that  the 
energy  will  enter  the  ocean  normal  to  the  surface.  The  field  at  some  point  in  the 
ocean  will  thus  be  approximately  of  the  form 


-ik  (y+d)  +ik^(y+d) 


+  p,  e 


(6) 


where 


UJ 
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(7) 


The  subscript  t  indicates  that  this  is  the  magnetic  field  in  the  ocean  due  to 
the  energy  entering  the  top  surface  only.  The  quantity  is  the  reflection  coeffi¬ 

cient  of  the  magnetic  field  at  the  bottom  of  the  ocean  and  is  given  by 


For  the  values  of  the  parameters  involved  here  is  very  nearly  equal  to  -  1.  Ih  (6) 
the  coefficient  A  is  chosen  so  that  the  magnetic  field  is  continuous  across  the  air-sea 


A. 


boundary.  Thus 


ik  d  +ik.d 

pbe 


(9)  > 


Contribution  Through  Land-Sea  Boundary. 

In  this  section  we  find  the  magnetic  field  in  the  ocean  due  to  the  energy  entering 
the  Land-Sea  boundary  only. 

To  first  obtain  an  approximate  expression  for  the  magnetic  field  at  the  bottom  sur¬ 
face  of  the  ocean,  due  to  the  energy  entering  this  interface,  Green’s  Theorem  will  be 
used. 

Consider  the  region  bounded  by  the  closed  dotted  line  path  1*  in  figure  2.  It  can 
2 

be  shown  that  the  magnetic  field  at  some  point  x,y  in  this  region  is  given  by 

h(r)  '  }  ^  [«rU'>  h(r')  -h(r->  ^  G(r|r*)  ]  dl'  CIO) 

where  G(r|r’),  the  Green's  function,  is  a  solution  to  the  inhomogeneous  wave  equation 


(V2  +  k2  )  G(r | r*)  «■  iu>e  6(|r  -  r*|) 


k2  -  )  «  y  cr2  (12) 

This  Green’s  function  represents  the  magnetic  field  at  r  due  to  a  line  source  at  r’ 
described  by  the  Dirac  delta  function  6(|r  -  r*|).  The  operator  -jj—  denotes  the 
derivative  normal  to  the  boundary  1*. 

The  Green’s  function  will  be  chosen  so  that  its  normal  derivative  will  vanish 

along  y  -  -  d  for  all  x.  In  addition  if  we  assume  that  the  ocean  is  a  perfect  con- 
3h 

ductor,  — ,  vanishes  for  y  =  -  d,  x  >  0.  In  this  case  the  integral  along  y  »  -  d, 
0  £  x  s  co  will  be  zero.  By  introducing  small  losses  the  .  -ral  over  the  infinite;. 


arc  will  also  be  zero. 


5. 


Thus  (10)  reduces  to 
The  solution  to  (11) 
that  ~  G(r|r)  ■  0  at  y  ■ 


an  integral  along  the  y  axis  from  y  «  -  «  tc  y  ■  -  d. 

for  G(r|r»)  along  y  -  -  d  subject  to  the  condition 

(3) 

-  d  is  well  known  and  is  given  by 


G(r| r') 


!IiH(i) 
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J (x-x*)2  +  (y,+d)2  ) 


\13) 


In  order  to  evaluate  (10)  we  now  need  to  know  the  magnetic  field  h-along  the  path 
of  integration.  This,  of  course,  is  not  known  exactly.  As  is  often  the  case  in  prob¬ 
lems  of  this  type  we  shall  approximate  h  over  this  path  by  assuming  it  is  equal  tc  the 
incident  magnetic  field  given  in  (3). 

Substitution  of  (13)  and  (3)  into  (10)  yields  the  following  expression  for  the 
magnetic  field  along  the  bottom  surface  of  the  ocean 


This  integral  cannot  be  readily  evaluated  for  all  values  of  x  However,  for 

j  2 

sufficiently  large  x,  say  for  x  >  6  6^  where  m  J  ^  ~  is  the  skin  depth  in  the 

o  2 

earth,  we  show  in  the  appendix  that  (14)  is  approximately  equal  to 
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(15) 


We  now  drop  our  assumption  of  infinite  conductivity  anr  write  that  the  magnetic 
field  in  the  sea  is  approximately  of  the  form 


i 
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6. 
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(16) 


The  subscript  b  indicates  that  this  is  the  magnetic  field  in  the  ocean  due  co 
the  energy  entering  through  the  bottom  of  the  ocean  only. 

The  reflection  coefficient  at  the  top  of  the  ocean  is  P,.  . 


.  ±1l2-U!L 


(17) 


+  -a 


J  ei  "  iUJ 

This  reflection  coefficient  is  also  nearly  equal  to  -1  for  our  purposes. 
The  coefficient  B  is  given  by 


B  .  _ _bs _ 

'  -ik  d  +ik  d 

e  +  pt  e 

which  makes  the  magnetic  field  continuous  at  the  bottom  of  the  sea. 


(18) 


Total  Magnetic  Field  in  the  Ocean 

The  total  magnetic  field  in  the  ocean  la  given  by  the  sum  of  (6)  and  (16). 

Substituting  for  A  and  B  from  (9)  and  (18)  we  obtain 

-ik^y+d)  +ik1(y+d)  ik.y  -k,v 

e  +  PK e  e  +  p„e  1 

L  ^  V  ^  — 1^ L  i 

tot.  „  -ik.d  +  ik.d  ts  -ik  d  +ik,d  *bs  1 

e  *  be  e  +  pfce 

If  we  set  p  and  equal  to  -1  we  have  finally 

sin  k^(y  +  d)  sin  k^y 

\otal  ~  sin  k^d  ^ts  “  sin  k^d  Njs 

where  h  and  h^s  are  given  by  (5)  and  (15)  respectively. 
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Electric  Field. 


Now  that  we  have  an  expression  for  h,  the  electric  field  is  easily  obtained 
from  (1).  In  particular  the  horizontal  component  of  the  electric  field  is 
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COS  k.  (y  +  d) 

— - -  h 

sin  k^d  ts 


cos  k,  y 

— — ~r  h, 

sin  k^  d  bs 


(21) 


where  again  we  are  making  use  of  the  fact  that  «  1. 

Numerical  Results. 

Numerical  values  for  the  x  component  of  the  electric  field  were  determined  by 
digital  evaluation  of  equation  (21),  The  results  are  plotted  in  Figs,  3  to  12  as  a 
function  of  depth  and  distance  from  shore  tor  frequencies  of  10,  100,  ar,d  1000  c/s. 
Depths  range  from  50  meters  to  800  meters  and  distances  from  shore  range  from  3  km 
to  40  km  .  The  depth  is  plotted  as  the  abcissa  and  the  various  values  of  x 
(distance  from  shore)  determine  a  family  of  curves. 

The  contribution  to  the  field  due  to  the  energy  passing  through  the  air-sea 
boundary  only  is  shown  as  a  dotted  line.  This  enables  one  to  see  more  explicitly  the 
effect  that  the  energy  passing  through  the  horizontal  land-sea  boundary  has  on  the 
total  field. 

In  Figs.  13  and  14  6^  and  $2>  the  skin  depth  in  the  sea  and  land  respectively 
are  plotted  as  a  function  of  frequency.  This  allows  interpretation  of  the  curves  in 
terms  of  skin  depth  instead  of  absolute  distance. 


Discussion. 

Before  discussing  the  results  it  must  be  admitted  that  both  the  model  and  the 
incident  field  which  we  have  chosen,  and  the  manner  in  which  they  approximate  the 
actual  situation  leave  much  to  be  desired.  Nevertheless,  the  approximations  used 
were  felt  to  be  sufficient  to  obtain  a  first  order  solution  to  our  problem,  namely  to 
establish  the  relative  importance  of  the  contributions  from  the  top  and  bottom  of  the 


ocean. 


Upon  examination  of  the  curves,  it  can  be  seen  that  the  contribytion  to. the 

net  field  which  enters  through  the  ocean-land  interface  is  quite  prominent  in  the 

region  near  the  shoreline.  As  the  distance  from  the  shore  increases,  this  contri^ 

bution  becomes  less  and  less  noticeable.  This  is  due  to  the  fact  that  k  ,  the  wave 

x 

number  which  governs  the  propagation  in  the  x  direction  for  the  contribution  through 
the  air- ocean  interface,  has  a  small  imaginary  part,  whereas  the  attenuation  of  the 
contribution  through  the  ocean-land  interface  is  quite  high  due  to  the-  conductivity 
of  the  earth. 

The  effect  of  the  depth  of  the  ocean  on  the  electric  field  is  also  quite  notice¬ 
able.  For  the  region  near  the  shoreline  the  curves  show  the  earth  propagated  contri¬ 
bution  more  distinctly  as  the  depth  of  the  ocean  increases.  This  is  understandable 
since  for  a  deep  ocean  the  energy  entering  the  top  surface  will  be  greatly  attenuated 
before  reaching  the  bottom.  The  contribution  through  the  bottom  will,  therefore,  be 
more  evident. 

Since  an  approximation  was  used  to  simplify  the  arguments  of  the  Hankel  function 
(see  appendix)  our  expressions  are  not  valid  for  x  <  6  6^.  The  model  itself  is,  however, 
not  very  reasonable  in  this  region. 


Appendix  A  • 

h^g  for  large  X  *r-r  ■ 

Consider  equation  (14) „  Because  of  the  decaying  exponential  character  of  the  inte¬ 
grand  ^2y  is  comP^ex)  the  contribution  to  the  integral  is  appreciable  only  for  that 
portion  of  the  integration  path  within  a  few  skin  depths  of  the  bottom  of  the  ocean.  The 
skin  depth  in  the  earth  is  given  by 


./  “V2 


(Al) 


If  we  assume  the  integrand  to  be  negligible  for  |y'[  greater  than  say  26^,  and 
restrict  ourselves  to  values  of  x  greater  than  66 ^  then  the  argument  of  the  Hankel  func¬ 
tion  can  be  greatly  simplified  over  that  portion  of  the  integrand  that  is  most  important; 

V 

i  .e. 


i  o  c 

J  x  +  (y*  +  d)  m  x 


(A?) 


Since  the  values  of  d  considered  are  much  less  than  b^  its  presence  in  the  argument 
does  not  affect  the  approximation  to  any  great  extent. 

With  this  approximation  the  Hankel  functions  are  no  longer  functions  of  y *  and  may 
be  moved  outside  of  the  integral  sign  leaving  a  simple  exponential  term  to  integrate. 
Upon  integrating  we  obtain 

ik2yd 

h.  «  -Hf -  k  H^1}  (k~x)  +  ik,  H(l)(k,x)  (A3) 

bs  2k„  x  o  2  2  1  2 
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electromagnetic  energy  propagation  via  air  into  the  ocean  versus  propagation 
through  the  soil  and  ocean  bottom  into  the  sea.  Vertical  profiles  ot  the 
horizontal  electric  field  are  plotted  for  several  values  of  frequency,  ocean 
depth,  and  the  distance  from  the  shoreline.  The  relative  importance  cf  the  soil 
path  is  seen  to  be  greater  for  shorter  distances  from  the  shore  and  for  deeper 
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